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Approximate  Confidence  Intervals  for  an  Exponential 
Parameter  from  a Sequential  Life  Test 

Abstract 

-i! 

A sequential  life  test  for  the  exponential  location 
parameter  was  given  by  Epstein  and  Sobel  (1955) . This 
sequential  test  may  be  modified  by  truncating  the  test  at 
rg  failures  and/or  at  total  test  time  t'  . There  may  be 
a need  or  a desire  to  also  estimate  the  parameter  after  the 
test  decision,  using  the  test  data.  Bryant  and  Schmee  (1979) 
have  given  confidence  intervals  for  the  mean  lifetime,  (jSjT, 
from  a truncated  sequential  test  scheme,  using  methods  which 
depend  heavily  on  numerical  techniques  using  a computer.  A 
more  flexible  approach  is  considered  using  a martingale 
inequality  which  was  also  given  by  Wald  (1947)  in  another 
context.  Interval  estimates  are  found  which  are  functions  of 
a positive  constant  d which  must  be  chosen  less  than  an  upper 
bound  which  is  itself  a function  of  the  number  of  failures 
observed.  It  is  suggested  that  d be  chosen  as  a function 
of  the  sample  path  (i.e.,  after  the  test  is  complete).  The 
validity  of  the  confidence  coefficient  comes  into  question  if 
this  posterior  selection  of  d is  employed.  Simulation  studies 
indicate  that  the  resulting  intervals  are  usually  conservative. 
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Approximate  Confidence  Intervals  for  an  Exponential 
Parameter  from  a Sequential  Life  Test 

I . Introduction 

The  theory  of  sequential  tests  was  developed  by 
Abraham  Wald  (1947)  . Epstein  and  Sobel  (1955)  applied  this 
theory  directly  to  get  a sequential  life  test  in  the 
exponential  case.  This  requires  the  assumption  that,  if 
t is  the  lifetime  under  consideration,  the  probability 
density  function  of  t is  given  by 

r 


f(T)  = 


-T/e 


e e 


if  t > 0 


^0  elsewhere, 

where  0 > 0 is  the  average  lifetime  of  an  item.  In  the 
sequential  exponential  life  test,  the  continuation  region  is 
bounded  by  two  parallel  lines  inthe  failure-total  time  1 ane 
as  depicted  in  figure  1.1.  Initially  n items  are  placed  on 
test  with  or  without  replacement.  At  each  failure  time  a 
decision  is  made  to  accept  or  reject  the  null  hypothesis,  or 
to  continue  the  test.  The  test  continues  as  long  as  the  sample 


Figure  1.1.  Sequential  life  test  boundaries  with  time  and 
failure  truncation. 


path  remains  inside  the  continuation  boundaries.  A 
modification  of  this  test  is  to  truncate  the  continuation 
region  at  r^  failures  and/or  at  total  time  t'  . This 
truncation  is  also  illustrated  in  figure  1.1.  These 
truncated  sequential  tests  are  used  by  the  Navy  as  reliability 
acceptance  tests  in  MIL-STD-781C  (1977) . The  tables  given  in 
MIL- STD- 7 8 1C  are  for  testing  0 = 0^  against  0 = 0^ 

where  0^  < 0^  . Since  values  of  0 considerably  different 
from  0q  or  0^  could  lead  to  the  acceptance  or  rejection  of 
Hq  , it  may  be  of  interest  to  also  estimate  0 following  the 
test  decision  using  the  test  data.  This  report  gives  a method 
for  finding  approximate  confidence  intervals  for  0 from 
the  sequential  test  or  the  truncated  sequential  test. 

Bryant  and  Schmee  (1979)  have  found  confidence  intervals 
for  0 following  the  truncated  sequential  tests  in 
MIL- STD- 7 8 1C . Their  results  require  the  application  of  a 
technique  developed  by  Aroian  (1963),  and  Epstein,  Patterson 
and  Qualls  (1963).  This  technique,  which  Aroian  calls  the 
"direct  method,"  reduces  the  continuation  region  to  a grid 
of  discrete  points,  including  a set  of  points  on  the  boundary 
of  the  continuation  region.  Using  an  iterative  procedure  and 
the  properties  of  the  exponential  distribution,  a probability 
for  each  point  on  the  boundary  can  be  found  as  a function  of 
0 . Confidence  intervals  for  0 can  then  be  found  using 
numerical  techniques  with  the  computer.  A drawback  of  this 
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method  is  the  amount  of  computation  required  and  its 
dependence  on  the  location  of  the  continuation  boundary. 

All  computations  may  have  to  be  repeated  if  the  value  of 
a test  parameter  is  changed. 

The  objective  of  this  report  is  to  find  a method 
which  is  easier  to  apply  computationally  and  perhaps  is 
more  versatile.  Our  approach  is  to  use  density  ratios 
defined  at  each  failure  time  and  apply  an  inequality  which 
was  given  by  Wald  (1947)  . He  suggested  the  inequality 
could  be  applied  to  sequential  estimation  problems.  O'Brien 
(1973)  uses  this  inequality  to  find  sequential  confidence 
intervals  for  the  shape  parameter  of  the  gamma  distribution. 

The  inequality  is  also  known  to  be  a special  case  of  a well 
known  martingale  inequality.  Robbins  (1970)  and  Lai  (1976) 
also  have  results  in  sequential  estimation  from  this  inequality. 

In  this  report,  the  inequality  is  applied  in  the  manner 
used  by  O'Brien,  wherein  the  density  ratio  involves  a positive 
constant  which  must  be  chosen  prior  to  the  experiment  in  some 
optimal  fashion.  The  derivation  of  the  intervals  in  this  report 
produces  a constraint  on  this  constant  which  is  a function  of 
the  number  of  failures  observed.  Consequently,  it  is  suggested 
that  the  intervals  may  be  improved  if  the  constant  is  chosen 
as  a function  of  the  sample  path  generated  by  the  failure  times. 
Whether  this  posterior  method  of  selecting  the  constant 
invalidates  the  desired  confidence  coefficient  remains  an  open 
question.  Computer  simulations  indicate  that  the  resulting 
intervals  remain  conservative,  and  that,  while  their  width  is 


-5- 


generally  greater  than  the  Bryant-Schmee  intervals,  the 
comparison  is  favorable  in  light  of  the  utility  of  these 
intervals . 


II . Testing  without  Replacement 

If  n items  with  exponential  lifetimes  are  placed 

on  test  without  replacement,  then  denote  the  ordered 

failure  times  as  t.  , t„ , ...»  t . Define  X.  s t.  - t.  , , 

3 2 n 1 i l-l 

where  t^  = 0.  The  , i = 1,  2,  . ..,  n,  are  independent 
with  density  functions. 


f (x^e) 


r 

(n-i+1) 

6 

0 


e 


- (n-i+1) x . 

0 1 

/ 


x.  > 0 
i 


elsewhere. 


Consequently  we  can  write  the  joint  density  of  X^,  X^,  ... 

X,  as 
k 


k 

n f ( x . ; 0 ) 
i=l  1 


or  as 


f (xlfx2, . . . ,xR) 


Z (n-i+l)x. 
n!  9 i=l  1 

k e 

(n-k) 16 


x.  > 0 

i 

i = l 1 2 , 


,k 


0 


elsewhere 


Equivalently 


fU^  x2. 


-Tk 


xk)  = 


n! 


(n-k)  ! 0 


where 


k k 

T,  = £ (n-i+l)x.  = £ t.  + (n-k)  t, 

k i=l  1 i=l  1 k 

= total  time  on  test  at  t,  . 

k 


Also,  define 


Y. 

l 


'£  (X  ,X2,.  . 
f (XX,X 


,x. ;d19) 


,X. ; 6 ) 

l 


if  i 


1,  2, 


Y 

n 


if  i > n. 


n 


At  this  point  d^  is  an  arbitrary  positive  constant.  It 
can  be  shown  quite  easily  that  Y^ , Y2 , ...  is  a martingale. 


E (Y 


k+1 


|Yk,...Y1]  = 


[k+1 

iV'W1 


k+1 

n f (X. ;6) 
i=l  1 


>Yk Y1 


k 

n f(X.;d.9) 
i=l  1 1 

k 


f(Xk+l;dl6) 

f(Xk+l;0) 


I 
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V 


rf,iwdi°>i 

oo 

f 

— v 

_ £(xk+li8) _ 

k 

0 

£(xk+l;dl6) 

£(xk+l,(,) 


f (xk+l'e) dxk+l 


Thus,  Y^ , Y^ , ...  is  a martingale  by  definition.  We  can 
now  use  a well  known  martingale  inequality  (see  Doob  [1953]) 
which  states  that,  for  e >0, 

(v 

P [max  Y.  > e ] E , for  all  k. 

i<k  1 e 

Since  E(Y^)  = 1,  the  bound  does  not  depend  on  k and  the 
inequality  becomes 


P [max  Y^  > e ] 


t 


or 

P[Yi  < e,  for  all  i]  > 1 - ^ . (2.1) 

We  note  that  the  inequality  ( 2.1  ) was  also  given  by 
Wald  (1948),  without  relying  on  the  martingale  property, 
in  the  following  lemma. 

Wald  Lemma:  Let  X^  be  a sequence  of 

random  variables  with  joint  density  f^(x^,...,x  ),  n = 

1,  2,  3,  ...,  under  hypothesis  H.  and  f.  (x. ,...,x  ) under 

i u n x n 

hypothesis  . Let  a be  a constant  between  0 and  1.  Then, 
under  HQ, 
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test 


Pf(fln//fOn)  < a ' for  a11  n]  2 1 “ a> 


In  the  context  of  the  sequential  exponential  life 


Y . 
1 


6ke  dl6 


(dj^O)  e 


Tk  — 

-k  I 'V1' 

dl  e 


Then,  applying  inequality  (2.1  ) , we  get 


Tk  ± 

-k  --T  <dr1J  i 

P[dL  K e 1 < i , for  all  k]  > 1 - a. 


Equivalently , 


Tk{1'dl) 

> -klnd^  + In  a,  for  all  k 


> 1 - a , 


The  right  side  of  the  inequality  in  brackets  is  positive 
if  d^  < Finally  we  have,  for  d^  < 


e < 


Tk(1'di) 


d^  In 


(a 


for  all  k 


> 1 - a. 


<2.2  ) 
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We  will  refer  to  this  upper  bound  on  0 as  U^. 
Inequality  (2.2  ) can  aloo  be  expressed  as 

P(G  < min  U,  ) > 1 - a. 

k K ~ 


Likewise  define 


Z . 
l 


f ( , X2 ' * * 

f (x1,x27. 


. ,X.;0/d2) 

. . ,xi; 0) 


if  i 


1,2 


, n 


Z 

n 


if  i > n, 


where  d2  is  an  arbitrary  positive  constant.  From  inequality 
( 2.1  ) we  have 

P [ Z,  < — , for  all  k]  > 1 - a, 
k a — 


where 


Zk  = 


d-,0] 


k -TRd2/0 


-Tk 

,k 

= d2  e 


This  yields 


k -0^(d2_1)  1 

P[d*  e < — , for  all  k] 


1 - ct, 


or 
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Tk(1_d2)  k 

P[ 5 < - In  (ad~),  for  all  k]  > 1 - a. 

w z — 


k 1/k 

If  - In  (ad2)  > 0,  i.e.,  d2  < 1/a  , then 


Tk<l-d2) 

-In  (adk ) 


0,  for  all  k 


> 1 - a. 


( 2.3  ) 


We  shall  refer  to  this  lower  bound  on  9 as  L^.  If  is  to 

1/k 

be  positive  it  must  be  true  that  d2  < 1 and  d2  < 1/a  ' . 

1/k 

Since  a < 1 , we  have  the  single  constraint  d2  < 1.  Here 
also,  it  is  true  that 

P [max  L.  < 8]  > 1 - a, 
k K 


since  (2.3  ) is  valid  for  all  k.  Applying  the  Bonferroni 
inequality  gives 

P (max  L,  <8 , min  U,>0)  > P (max  L,<0)  + P (min  U.  >0)  - 1, 
k K k K k K k K 


or 


P [max  L.  < 9 < min  U,  ] >l-a  + l-  a-  l = l-2a. 
k K k K 

This  establishes  the  following  theorem. 

Theorem  2 . 1 If  a sequential  exponential  life 
test  is  conducted  without  replacement  and  terminates  with 
r failures,  then  a two-sidea  confidence  interval  for  the 


-1.1  - 


mean  lifetime  0,  having  confidence  coefficient  at  least 
1 - 2a,  is  given  by 


T (l-d2)  Tk(l-d  ) 

max  t— - min  

k<^r  -ln(ad2)  ' k£r  d^ln 

dl 


1 /k 

where  0 < d^  < a ' and  0 < d^  < 1. 

Due  to  the  nature  of  the  interval,  it  may  happen 

that 


max  min  . 

k<r  k - k<r  uk* 

This  leaves  us  with  a confidence  interval  whose  lower  limit 
exceeds  its  upper  limit.  Computer  simulation  studies  have 
shown  that  this  inversion  of  the  endpoints  occurs  infre- 
quently, often  as  a consequence  of  an  unusually  early 
failure.  To  avoid  this  predicament,  the  following  iterative 
scheme  is  proposed  whenever 

max  T min  .. 
k<r  ^k  - k<r  k* 

(1)  The  upper  limit  of  the  confidence  interval,  say 
U,  is  assigned  the  value  of  the  next  larger 
member  in  the  set  {U^  : k = l,2,...,n}. 


(2)  If  U is  greater  than  the  lower  limit,  say  L, 
the  endpoints  are  ordered  correctly  and  the 
procedure  stops. 

(3)  If  U is  smaller  than  the  lower  limit,  assign  to 
L the  value  of  the  next  smaller  number  of  the 
set  {L^  : k = 1,2,. ..,n}. 

(4)  If  L > u return  to  step  (1). 

Since  inequalities  ( 2.2  ) and  ( 2.3  ) are  valid  for  all  k, 
the  interval  that  results  from  applying  this  procedure  will 
still  have  a confidence  coefficient  at  least  l-2ot. 

Accept  decisions  in  the  sequential  test  occur 

between  failure  times.  Consequently,  if  an  accept  decision 

occurs  and  j failures  have  been  observed,  the  total  time 

on  test  when  the  decision  is  made  to  accept  Hq,  denoted  by 

T . say,  will  exceed  T.,  the  total  time  on  test  at  the  jth 
a»  J 3 

failure.  Furthermore 

-lnfad^)  -lntad^) 


and  it  may  also  be  true  that 


-In (ad^) 


2 


) 


max  Tk(1-d2) 


^ -lnfad^) 


This  suggests  that  we  might  shorten  the  confidence 
interval  for  0 by  increasing  the  lower  limit  utilizing  this 


total  test  time  at  the  time  of  an  accept  decision.  Note 

that  the  upper  limit  cannot  be  improved  in  the  same  manner 

since  there  we  seek  the  minimum  of  a set,  and  T . > T. 

a/1  3 

implies  the  minimum  will  occur  at  a failure  time.  The 


lower  limit 

of 

the 

interval 

in  Theorem 

2.1 

was  chosen  as 

the  maximum 

of 

the 

set  (L^, 

k = 1,2,... 

fr} . 

The  following 

lemma  admits  the  use  of 


Ta.i(1-d 

-In (ad^ ) 


as  the  lower  limit  when  an  accept  decision  occurs. 

Lemma  2 . 2 If  the  sequential  life  test  terminates 

in  an  accept  decision  with  j failures  and  total  test  time 

T then 
a / J 


P 


e > 


Ta,j(1-d2n 


-In (ad^) 


> 1 


a . 


Proof:  We  note  that  this  is  an  exponential  life  test, 

without  replacement,  which  is  observed  for  a total  test  time 

t'.  This  test  can  be  shown  to  be  equivalent  to  a Poisson 

process  with  parameter  X = 1/0  which  is  observed  for  actual 

time  t'.  We  will  apply  that  result  here  with  t'  = T ..  A 

a,  j 

well  known  result  for  a Poisson  process  (see  Barlow  and 
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Proschan  (1975),  Theorem  3.7)  may  bo.  adapted  as  follows: 

Given  j failures  in  a Poisson  process  in  (0,  T .], 

a / j 

the  distribution  of  the  failure  times  (here  interpreted  as 

actual  times)  , T2 , . ..,  T_.  is  the  distribution  of  the 

order  statistics  from  the  uniform  distribution  over  the 

interval  ( 0 , T . ] . 

a , j 

As  a consequence,  if  y^  represents  the  observed 
value  of  iv,  then 


f (y1/y2/ • • • /yj I j) 


i l 


a, 3 


/ 


and 


P(j  Failures)  = 


T < 

e 9 (T  ./0) 
_JLtJ 

-i  ! 


Thus 


f (y1*y2****  >Vy  j) 


Now  if  we  apply  Wald's  lemma  in  the  manner  leading  to 
Theorem  2.1  letting 
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y 


i 


u 


‘d2Ta, 

d2  e 


• /Q 


and  f 


-T  ./0 

a, 3 


0- 


Oj 


the  conclusion  of  the  lemma  follows  by  the  same  algebraic 
manipulations . 

There  remains  the  problem  of  selecting  the  values  of 

and  d2  which  are  in  some  sense  optimal,  or  near  optimal, 

while  satisfying  the  given  constraints.  We  shall  say  the 

values  of  d^  and  d2  are  optimal  if  they  minimize  the  length 

of  the  confidence  interval  in  Theorem  2.1  . The  constraint 
1/k 

d-^  < a ' , k = 1,  2,  ...,  r,  poses  a dilemma.  If  d^  is  a 

fixed  constant  as  the  theory  behind  Theorem  2.1  requires, 
then  it  must  be  true  that  d^  < a.  It  is  not  likely  that  a 
value  of  d^  in  this  range  will  be  "optimal"  in  any  sense  of 
the  word.  On  the  other  hand  if  a value  of  d^  is  chosen 
greater  than  a,  and  the  test  terminates  with  one  failure, 
then  there  is  no  upper  bound  (other  than  + °°)  for  the 
confidence  interval. 

The  martingale  inequality  on  which  Theorem  2.1 
is  based  has  wide  applicability.  It  is  reasonable  to 
suspect  that  in  a specific  application  the  bound  may  be 
quite  conservative.  Specifically  the  probability  in 
Theorem  2.1  may  considerably  exceed  l-2a.  This  suggests 
that  some  liberties  might  be  taken  in  selecting  d^  and  d2 
without  violating  the  probability  inequalities  in  ( 2.2) 
and  ( 2.3  )•  Therefore,  we  propose  tochoosc  d1  and  d2  as 
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functions  of  the  sample  path  generated  by  the  sequential 
test.  The  following  scheme  implements  this  proposal. 

Treat  U^,  k = 1,  2,  . ..,  n,  as  a function  of  and 
for  each  failure  time,  t^ , find  the  value  of  d^  which 
minimizes  U^.  This  optimal  value  is  found  to  be  the  root 
of  the  equation 

g(d)  = 1 - (i)  In  u -d  + lnd. 

Note  that  the  solution  is  not  dependent  on  T^  and  so,  for 

a given  a,  a set  of  optimal  values  of  d^,  say  fd^,  d^2* 

^13'  •**'  } Can  computed  by  numerical  methods.  A 

table  of  these  values,  for  several  values  of  a,  is  given  in 

Appendix  1.  At  each  failure  time,  t^,  compute  using  d^. 

Then  if  the  minimum  value  of  U.  occurs  at  t say,  we  set 

k mJ 

d^  E d^m<  The  following  lemma  shows  that  for  each  k,  the 

1/k 

relationship  d^k  < a ' will  hold. 

Lemma  2.3  . If  d^  is  the  root  of  the  equation 
g(d)  = 1 - (1/k)  In  a -d  + lnd,  then  dlk  < a1^. 

Proof:  First,  g’ (d)  = 1/d  - 1 > 0,  which  implies 

that  g(d)  is  increasing  on  the  interval  (0,  1).  Also 
g ( 0 ) = -°°  and  g(l)  =-(l/k)  In  a > 0.  It  follows  that  g(d) 
crosses  the  axis  once  from  below  in  (0,  1).  But  gfa'*'^^)  = 
l-a1/k  > 0.  Thus  dlk  < a1/k. 

In  the  same  way,  consider  as  a function  of  d2 
and  find  the  set  ^21'^22 ' * ' * f<^2n^  values  which  maximizes 
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for  k=lf  2,  . ..,  n.  These  values  are  the  roots  of  the 
equations  h(d)  = 1 - (1/k)  In  a - lnd  - d ^ and  are  also 
tabled  for  some  specific  values  of  cx  in  Appendix  1.  If  the 

| 

maximum  value  of  occurs  at  t^,  set  d2  H ^2?,’  Usin9  the  ! 

same  steps  as  in  lemma  2.3  , it  is  seen  that  d2k  < 1 for 

each  k.  I 

To  clarify  this  procedure,  consider  the  following  ; 

example:  Test  plan  I of  MIL-STD-781C  is  used  to  test 

Hq  : 0 = 0q  against  : 0 = 0^,  where  Og/O^  = 1.5, 
a = 3 = .1.  The  test  ends  with  an  accept  decision  based  on 
six  observed  failures.  At  each  failure  time,  the  total  time 
on  test  was  computed  as  were  the  upper  and  lower  95% 
confidence  limits  using  (2.2  ) and  (2.3  ) and  the  optimal 
values  of  d^  and  d2  given  in  Appendix  1 . These  computations 
are  summarized  in  Table  2.1  . A lower  confidence  limit  was 
also  computed  at  the  acceptance  boundary  as  described  by 
lemma  2.2,  using  T , = 13.91.  Failure  times  were 

a / O 

simulated  using  0 = 0q  = 1.5.  According  to  our  procedures 

based  on  Theorem  2.1  , an  approximate  90%  confidence 

interval  for  0 is  [.984,  4.101]  with  d^  = .2628  and  d2  = 

.4244.  The  Bryant-Schmee  90%  confidence  interval  for  this 

example  is  [1.108,  4.681]. 

The  second  method  is  really  a simplification  of 

this  first  method.  During  computer  simulations  it  was 

noticed  that,  for  a fixed  d.  or  d_,  max  L.  and  min  U, 

k<r  K k<r  k 
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TABLE  2 . I 

Computation  Results  for  Confidence  Interval  Example 


Failure 

k 

Total 

Test 

Time 

Tk 

dlk 

Uk 

d2k 

Lk 

1 

1.362 

.0187 

72.649 

.1741 

.2371 

2 

2.623 

.0900 

14.572 

.2601 

. 3411 

3 

4.618 

.1589 

9.690 

. 3181 

.4896 

4 

4.815 

.2159 

5.576 

.3615 

.4351 

5 

5.389 

.2628 

4.101 

.3959 

.4267 

6 

12.555 

. 3020 

6.928 

.4244 

. 8880 

a,  6 


A 


■ 13.91 


4244 


9838 
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usually  occurred  at  or  near  t .,  the  last  observed  failure 

/ * 

time.  Hence,  if  we  set  = d.  and  d~  = d-.  we  should  get 
results  nearly  the  same  as  with  the  first  method.  It  was 
also  felt  that  this  second  method  would  be  easier  to  justify 
analytically.  Unfortunately,  an  analytic  justification  of 
either  of  these  methods  remains  an  unsolved  problem.  It  is 
not  known  if  these  posterior  methods  of  selecting  d^  and  d£ 
invalidate  the  probability  inequalities  (2.2  ) and  (2.3  )• 

To  obtain  some  empirical  conclusions,  we  conducted 
computer  simulations  of  two  truncated  sequential  test  plans 
from  MIL-STD-781C.  At  the  time  of  test  decision  60%,  80%, 
and  90%  confidence  intervals  for  the  mean  failure  time  6 
were  formed  using  both  methods  discussed  here  for  selecting 
d^  and  d£.  These  confidence  coefficients  were  chosen  so 
that  our  results  could  be  compared  with  those  obtained  by 
Bryant  and  Schmee  (1979).  Each  simulation  used  500  trials. 
Estimates  of  the  confidence  coefficient,  and  the  average 
value  of  the  endpoints  were  computed  for  our  methods  and  for 
the  Bryant  and  Schmee  intervals.  Partial  results  are  given 
in  Tables  2.2  through  2.5  . It  was  noted  that,  when  the 
test  ended  with  an  accept  decision  and  one  failure,  the 
upper  endpoint  on  the  confidence  interval  was  often  very 
large,  which  tended  to  inflate  the  average  endpoints. 

Because  of  this,  the  average  endpoints  were  also  computed 
excluding  those  trials  which  terminated  in  an  accept 
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TABLE  2 . 2 


Simulation  Results  for  909.  Confidence 
MIL- STD- 7 8 1C , Test  Plan  I;  0/0 
a = 6 = .1  U 


Interval  from 
= 1.5 


0 

Estimated 

Coefficient 

Ave.  Endpoints, 
at  Least  1 Failure 

Ave.  Endpoints, 
at.  Least  2 Failures 

. 884 

. 307 

r 

. 917 

. 307 

9 

.917 

.5 

. 898 

.259 

/ 

1.059 

.259 

9 

1.059 

. 946 

.256 

9 

1.102 

.256 

9 

1.102 

. 912 

.714 

9 

2.399 

.708 

9 

1.637 

1 

. 917 

.655 

9 

3.208 

.650 

9 

1.760 

. 964 

.647 

9 

3.228 

.642 

9 

1.804 

. 916 

1.020 

9 

7.038 

1.003 

9 

3.743 

1.5 

.912 

.942 

9 

8.487 

. 930 

9 

4.192 

. 962 

.941 

9 

8.553 

.930 

9 

4.367 

.912 

1.137 

9 

15.331 

1.096 

9 

5.800 

2.0 

.870 

1.018 

9 

20.113 

.989 

9 

6.734 

.906 

1.017 

9 

20.265 

.988 

9 

7.081 

Note : 

In  each  of 

Tables 

2.2 

through 

2.5  , 

(1)  The  first  entry  in  each  cell  corresponds  to 
Bryant-Schmee  results. 

(2)  The  second  entry  corresponds  to  the  first 

method:  using  the  optimal  d.  and  d-  at  each 

failure . 

(3)  The  third  entry  corresponds  to  the  second 

method:  using  the  optimal  d.  and  d_  from  the 

final  failure  time. 


0>4*MSV< 
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TABLE  2.3 


Simulation  Results  for  90%  Confidence  Interval  from 
MIL- STD- 7 8 1C , Test  Plan  III;  Oq/0.  = 2, 
a = 8 = . 1 


e 

Estimated 

Coefficient 

Ave.  Endpoints, 
at  Least  1 Failure 

Ave.  Endpoints 
at  Least  2 Failures 

. 886 

.259 

9 

1.398 

.259 

9 

1.398 

.5 

.944 

. 195 

9 

1.733 

.195 

9 

1.733 

.974 

.193 

9 

1.803 

. 193 

9 

1.803 

.912 

.868 

9 

3.746 

. 558 

9 

2.396 

1 

.931 

.472 

9 

4.344 

.463 

9 

2.715 

. 962 

.469 

9 

4.401 

.460 

9 

2.797 

.938 

.792 

9 

8.729 

. 768 

9 

3.891 

1.5 

.944 

.707 

9 

11.239 

.688 

9 

4.272 

.970 

.706 

9 

11.430 

.687 

9 

4.521 

.956 

.940 

9 

18.898 

.895 

9 

5.838 

2.0 

. 962 

. 855 

9 

25.942 

.825 

9 

7.232 

.986 

.855 

9 

26.140 

.825 

9 

7.616 

TABLE  2.4 


Simulation  Results  for  80%  Confidence  Interval  from 
MIL- STD- 7 8 1C , Test  Plan  I 


0 

Estimated 

Coefficient 

Ave.  Endpoints, 
at  Least  1 Failure 

Ave.  Endpoints, 
at  Least  2 Failures 

.784 

. 342  , 

. 798 

. 342  , 

.798 

.5 

.816 

.286  , 

.907 

.286  , 

.907 

. 898 

.281  , 

.964 

.281  , 

. 964 

. 810 

.760  , 

1.871 

.752  , 

1.504 

1 

.845 

.696  , 

2.434 

.691  , 

1.611 

. 920 

.680  , 

2.492 

.674  , 

1.696 

.808 

1.125  , 

4.485 

1.103  , 

3.022 

1.  5 

. 790 

1.018  , 

5.323 

1.003  , 

3.456 

.892 

1.011  , 

5.385 

.996  , 

3.609 

.867 

1.290  , 

8.433 

1.247  , 

4.661 

2.0 

.802 

1.120  , 

10.425 

1.091  , 

5.298 

.867 

1.118  , 

10.661 

1.089  , 

5.657 

*r 
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TABLE  2.5 

Simulation  Results  from  801  Confidence  Interval  from 
MIL-STD-781C,  Test  Plan  III 


0 

Estimated 
Coef  f ici ent 

Ave.  Endpoints, 
at  Least  1 Failure 

Ave.  Endpoints 
at  Least  2 failures 

. 788 

. 294  , 

1.085 

. 294 

9 

1.085 

.5 

.872 

.215  , 

1.346 

.215 

9 

1.346 

.922 

.212  , 

1.4  37 

.212 

9 

1.437 

. 830 

.601  , 

1.857 

.601 

9 

1.857 

1.0 

. 887 

. 489  , 

2.098 

.489 

9 

2.098 

.920 

. 484  , 

2.198 

.484 

9 

2.198 

. 856 

.916  , 

5.356 

.882 

9 

2.978 

1.5 

.849 

.798  , 

7.064 

.772 

9 

3.527 

.900 

.794  , 

7.188 

.769 

9 

3.731 

.892 

1.081  , 

9.882 

1.020 

9 

4.311 

2.0 

.912 

. 946  , 

12.656 

.903 

9 

5.043 

.956 

.945  , 

12.693 

.903 

9 

5.332 

-24- 


decision  with  one  failure.  As  expected,  the  Bryant-Schmee 
intervals,  on  the  average,  were  somewhat  shorter  than  the 
intervals  from  both  of  our  methods.  However,  we  also  noted 
that  for  some  sample  paths,  our  methods  gave  shorter 
intervals  than  the  corresponding  Bryant-Schmee  intervals. 

The  simulation  results  also  indicate,  that  the  intervals 
generated  from  our  posterior  methods  of  selecting  d^  and  d^ 
remain  conservative.  We  also  note  that  our  methods  may  be 
just  as  easily  applied  to  an  untruncated  sequential  test, 
while  this  is  not  so  with  the  Bryant-Schmee  approach.  Their 
method  requires  non-trivial  computation  of  probabilities  at 
every  discrete  acceptance  point  up  to  that  point  where  the 
sample  path  ended.  In  the  untruncated  sequential  test  the 
sample  path  may  continue  for  a long  time,  resulting  in  a 
rather  formidable  computation  problem  to  get  the  Bryant- 
Schmee  intervals.  We  believe  that,  considering  the 
relative  ease  of  computation  and  application,  these  approxi- 
mate confidence  intervals  provide  a competitive  alternative 
to  those  of  Bryant  and  Schmee. 

Ill-  A Confidence  Interval  with  No  Observed  Failures 

The  sequential  life  test  may  end  in  acceptance  with 
no  failures  observed.  In  this  situation  the  confidence 
interval  should  be  one  sided  since  there  is  no  information 
on  which  to  base  an  upper  limit.  If  no  failures  are 
observed  by  total  time  Tq,  then  we  have  observed  a Poisson 
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process  with  parameter  X = n/9  for  a total  time  Tg , or 
actual  time  Tg/n.  This  is  equivalent  to  type  I,  or  time, 
censoring  and  Epstein  (I960)  showed  that  a one-sided  100 (1-a) 
percent  confidence  interval  for  0 is 


J 


(3.1) 


When  no  failures  are  observed,  the  two-sided 
confidence  interval  of  Theorem  2.1  reduces  to 


T0(l-d2)  T0(l  V 


-In  a ' dj  In  a 


( 3.2  ) 


If  k = 0,  the  constraint  on  d^  becomes  0 < d^  < 0,  leaving 
the  upper  endpoint  undefined.  Since  no  failure  times  are 
known,  we  define  the  upper  endpoint  to  be  infinity.  The 
lower  endpoint  decreases  as  d ^ decreases  and  approaches 
Tg/(-ln  a)  as  goes  to  zero  in  the  limit.  The  resulting 
interval  is 


- In  a 


(3.3  ) 


But  from  the  relationship  between  the  chi-square  and  Poisson 


distributions  we  have 


- 
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Hence , 

~ ln  a = Kf2/2  ' 

and  ( 3.1  ) is  equivalent  to  ( 3.3  ),  and  either  provides  a 
one-sided  100  (1-a)  percent  confidence  interval  for  0 when 
no  failures  have  been  observed  in  the  life  test. 

1 V . Testing  with  Replacement 

In  this  section  we  will  show  that  the  derivations  in 

Section  TI  are  also  valid  if  the  sequential  test  is 

conducted  with  replacement.  Testing  with  replacement  is 

equivalent  to  observing  a Poisson  process  with  parameter 

A = n/0 , where  again,  n is  the  number  of  items  placed  on 

test.  As  in  Section  II  , define  X.  = t.  - t.  . . The  X., 

l l l-l  l 

i = 1,  2,  ...»  n,  are  independent  and  are  exponentially- 
distributed  with  mean  lifetime  0/n.  Hence, 


f (x^ 


• / Xj^ ; 0 ) 


where  T^  = nt^  = total  time  on  test  at  t^.  We  see  that  the 
joint  densities  of  X^,  X^ , ...,  X^  when  testing  with  or 
without  replacement  differ  only  by  a constant  involving  n 
and  k.  Thus,  the  density  ratios  in  Section  II  will  be  the 
same  as  here  in  the  replacement  case.  It  follows  that 


Theorem  2.1  also  holds  when  testiny  with  replacement. 

Now,  to  establish  the  analogous  result  to  Lemma  2.2  , 

we  note  that  a total  test  time  of  T is  equal  to  actual 

a,  D 

time  t . = T ./n  if  failed  items  are  replaced.  If  t.  is 
a, j a, j c i 

the  ifc^  ordered  failure  time  and  y^  represents  the  observed 
value  of  t . , then 

l 

f<yi'V2 Vj 1 3)=  - — j * 

ta<  j 


and 


nt 


AsJ. 


P(j  Failures)  = 


nt 

_a,_2 

e 


j! 


Thus 


f (y1*y2 


nt 

j SLU 

nJ  0 

'3)  = — r e 

e3 


T 

-_a-q 


n- 


Except  for  the  constant  n3 , this  is  the  same  result 
obtained  in  the  proof  of  Lemma  2.2  . The  conclusion  of 

Lemma  2.2  , for  testing  with  replacement,  follows  in  the 
same  manner  as  in  that  proof. 

It  follows  that  all  results  and  conclusions  of 
Section  II  , where  testing  is  without  replacement,  also 
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apply  to  the  sequential  test  scheme  where  failed  items  are 
replaced . 
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APPENDIX  1 

Tabled  here  are  the  optimal  d 2 values  for  determining 
the  lower  limit  of  the  100 (1-a)  percent  two-sided  confidence 
interval . 


Failures 

k 


Confidence  Coefficient 
99i.  9 8%  95%  90%  80% 


1 

0.1186 

0.1309 

0.1522 

0.1741 

0.2045 

2 

0.1879 

0.2045 

0.2323 

0.2602 

0. 2973 

3 

0.2380 

0.2568 

0.2876 

0. 3180 

0. 3575 

4 

0.2771 

0.2973 

0.3298 

0. 3615 

0.4023 

5 

0. 3092 

0.3301 

0.3637 

0. 3962 

0.4371 

6 

0. 3364 

0.3575 

0. 3918 

0.4243 

0.4656 

7 

0.3597 

0.3812 

0.4160 

0.4485 

0.4894 

8 

0. 3804 

0 . 4023 

0.4366 

0.4691 

0.5096 

9 

0.  3988 

0.4208 

0.4551 

0.4876 

0.5271 

10 

0.4151 

0.4371 

0.4718 

0.5034 

0.5430 

11 

0.4300 

0.4516 

0.4867 

0.5184 

0.5570 

12 

0.4437 

0.4656 

0.4999 

0.5315 

0.5702 

13 

0.4560 

0.4779 

0.5122 

0.5430 

0.5816 

14 

0.4674 

0.4894 

0.5236 

0.5544 

0.5922 

15 

0.4779 

0.4999 

0.5342 

0.4641 

0.6027 

16 

0.4885 

0.5096 

0. 5430 

0.5746 

0.6115 

17 

0.4973 

0.5192 

0.5518 

0.5825 

0 .6203 

18 

0.5061 

0.5271 

0.5605 

0.5904 

0.6273 

19 

0.5148 

0.5359 

0.5685 

0.5992 

0.6344 

20 

0.5219 

0. 54  30 

0.5764 

0.6062 

0.6414 

21 

0.5298 

0.5509 

0.5834 

0.6115 

0.6484 

22 

0.5368 

0.5570 

0.5904 

0.6186 

0.6537 

23 

0.5430 

0.5641 

0.5957 

0.6238 

0.6590 

24 

0.5500 

0. 5702 

0.6027 

0.6309 

0.6660 

25 

0.5553 

0.5764 

0.6080 

0.6  361 

0.6695 

26 

0.5614 

0.5816 

0.6133 

0.6414 

0.6748 

27 

0.5676 

0.5869 

0.6186 

0.6449 

0.6801 

28 

0.5729 

0.5922 

0.6238 

0.6520 

0.6836 

29 

0.5781 

0.5975 

0.6273 

0.6555 

0.6871 

30 

0.5816 

0.6027 

0.6326 

0.6590 

0.6924 

31 

0.5869 

0.6062 

0.6379 

0.6643 

0.6977 

32 

0.5922 

0.6115 

0.6414 

0.6678 

0.7012 

33 

0.5957 

0.6150 

0.6449 

0.6730 

0.7047 

34 

0.5992 

0.6203 

0.6484 

0.6766 

0.7082 

35 

0.6045 

0.6238 

0.6520 

0.6801 

0.7117 

— . , C . . 


•v 
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Failures  Confidence  Coefficient 


k 

99% 

o'.* 

CC 

95% 

90% 

80% 

36 

0.6080 

0.6273 

0 .6555 

0.6836 

0.7152 

37 

0.6115 

0.6309 

0.6590 

0.6871 

0.7152 

38 

0.6150 

0.6344 

0.6625 

0.6889 

0.7187 

39 

0.6186 

0.6379 

0.6660 

0.6906 

0.7223 

40 

0.6221 

0.6414 

0.6695 

0.6941 

0.7258 

41 

0.6256 

0.6449 

0.6730 

0.6977 

0. 7293 

4 


i 


* 
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Tabled  here  are  the  optimal  d^  values  for 
determining  the  upper  limit  of  the  100(1 -a)  percent  two- 
sided  confidence  interval. 


Failures  Confidence  Coefficient 


k 

99% 

9 8% 

95% 

90% 

80% 

1 

0.0018 

0.0037 

0.0093 

0.0187 

0.0382 

2 

0.0267 

0.0382 

0.0619 

0.0900 

0.1329 

3 

0.0673 

0.0864 

0.1215 

0.1590 

0.2108 

4 

0.1090 

0.1329 

0.1742 

0.2156 

0.2712 

5 

0.1478 

0.1742 

0.2190 

0.2625 

0.3190 

6 

0.1825 

0.2108 

0.2576 

0.3020 

0. 3590 

7 

0.2137 

0.2429 

0.2907 

0. 3356 

0. 391 2 

8 

0.2415 

0.2712 

0. 3190 

0. 36  39 

0.4195 

9 

0.2664 

0.2966 

0. 3444 

0.  389  3 

0.4439 

10 

0.2888 

0. 3190 

0. 3678 

0.4107 

0.4654 

11 

0. 3098 

0. 3405 

0.3873 

0.4312 

0.4849 

12 

0.3288 

0.3590 

0.4059 

0.4498 

0.5005 

13 

0.3464 

0.3756 

0.4225 

0.4654 

0.5161 

14 

0.3620 

0.3912 

0.4381 

0.4800 

0.5298 

15 

0.3766 

0.4068 

0.4517 

0.4927 

0. 5434 

16 

0. 3903 

0.4195 

0.4654 

0.5064 

0 . 5551 

17 

0.4029 

0.4322 

0.4771 

0.5181 

0.5668 

18 

0.4146 

0.4439 

0.4888 

0.5278 

0.5746 

19 

0.4264 

0.4556 

0.4985 

0.5376 

0.5844 

20 

0.4371 

0.4654 

0.5083 

0.5473 

0.5942 

21 

0.4458 

0.4751 

0.5181 

0.5571 

0.6020 

22 

0.4667 

0.4849 

0.5278 

0.5649 

0.6098 

23 

0.4654 

0.4927 

0.5356 

0.5727 

0.6176 

24 

0.4732 

0.5005 

0.5434 

0.5785 

0.6234 

25 

0.4810 

0.5083 

0.5512 

0.5864 

0.6293 

26 

0.4888 

0.5161 

0 . 557  j 

0.5942 

0.6371 

27 

0.4966 

0.52  39 

0.5629 

0.6000 

0.6410 

28 

0.5044 

0.5298 

0.5707 

0.6059 

0.6488 

29 

0.5109 

0.5376 

0.5766 

0.6117 

0.6527 

30 

0.5161 

0.5434 

0.5824 

0.6176 

0.6566 

31 

0.5239 

0.5493 

0.5883 

0.6215 

0.6605 

32 

0.5298 

0.5531 

0.5942 

0.6273 

0.6683 

33 

0.5356 

0.5610 

0.5981 

0.6332 

0.6722 

34 

0.5395 

0.5668 

0.6039 

0.6371 

0.6761 

35 

0.5454 

0.5707 

0.6098 

0.6410 

0.6800 

36 

0.5512 

0.5746 

0.6137 

0.6449 

0.6839 

37 

0.5551 

0.5805 

0.6176 

0.6488 

0.6878 

38 

0.5610 

0.5844 

0.6215 

0.6527 

0.6917 

39 

0.5649 

0.5903 

0.6254 

0.6566 

0.6956 

40 

0.5707 

0.5942 

0.6293 

0.6605 

0.6995 

41 

0.5746 

0.5981 

0.6332 

0.6644 

0.7034 
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